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Particle production from vacuum fluctuations during inflation is briefly revisited. The moduli problem oc- 
curring with light particles produced at the end of inflation is addressed, namely the fact that some results are 
in disagreement with nucleosynthesis constrains. A universal solution to this problem is found which leads to 
reasonable reheating temperatures in all cases. It invokes the assumption that, immediately after inflation, the 
moduli evolve like non-relativistic matter. The assumption is justified in the context of massive chaotic inflation 
were, at the end of inflation, the universe evolves as if it was matter-dominated. 
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1. Introduction. — Particle production from vacuum fluc- 
tuations during inflation has been considered for some time 
(see, for instance, HHH). However, for light particles the res- 
ults obtained have been reported not to be in agreement with 
constrains coming from nucleosynthesis. This is the so called 
moduli problem related to cosmological particle production. 
In this letter we will try to better understand this issue in an 
attempt to clarify the existing conflict between the results ob- 
tained in 031] and |0] concerning such moduli problem. Being 
more precise, for a light moduli field, in Ref. H an abnor- 
mally small reheating temperature was obtained while, quite 
on the contrary, in Ref. (01 an abnormally big one was repor- 
ted. 

We will here try to point out the basic differences between 
the approaches used in each of these papers and then propose 
a common solution to the moduli problem, in the context of 
massive chaotic inflation. To summarize the situation, it is 
our view that a common main assumption made in each one 
of the above references was to establish that, after inflation, 
the magnitude of the moduli field remains almost constant for 
a while, until it starts to oscillate. Assuming, on the contrary, 
that the moduli field actually evolves, after inflation, like non- 
relativistic matter, the problem is most naturally avoided in 
both cases. This is quite remarkable given the fact that the 
previous results where so extremely different in the two cases 
considered. To justify the new assumption, in the context of 
massive chaotic inflation one may argue that at the end of in- 
flation the universe evolves in fact like a matter-dominated 
one. This will be explained in detailed in what follows. 

Using the same notation as in Ref. [4], the massive inflaton 
field, with mass m, will be called <j>, and the moduli fields 
X, the corresponding mass being m x . We will also set a = 
rrv^/m 2 (notice that in 121 the mass of the inflaton is denoted 
by M and the one of the moduli, by m). 

2. The massive inflationary model. — In this model, the 



Friedmann equation is given by 
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where V = m 2 ^ 2 /2, M p ~ 10 18 GeV is the reduced Planck 
mass, and the inflaton field equation reads 
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In the slow-roll approximation, that is, when the following 
equations hold, 
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the evolution of the Hubble parameter is given by 
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Inflation ends here at time tt = 3Ho/m 2 . To obtain the ap- 
proximate value of the Hubble parameter at the end of infla- 
tion, we impose that in this epoch the slow-roll parameters, 
namely 
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be of order one. Then, one has <j>(tf) ~ M p , which implies 
that H(tf) ~ Tin. Once the inflaton field strength drops below 
the corresponding Planck mass value, it begins to oscillate and 
the universe expands, like in the case of a matter-dominated 
universe (see, e.g., |5[], pgs. 235-43). 

On the other hand, assuming that the moduli field is minim- 
ally coupled to gravity, then the equation for the inhomogen- 
eous Fourier modes reads 
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Xk = 0. 
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Introducing now the variable \ — a X> ar, d using conformal 
time, this equation can be written as 
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The energy density, p x , is given by 
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This quantity is divergent and needs to be re-normalized. One 
way to do that is to use adiabatic substraction, which can- 
cels the ultraviolet modes and thus in fact only infrared modes 
contribute to the energy density. Actually, we are only inter- 
ested in modes which leave the horizon during inflation. Thus, 
at the end of the inflation era we only consider modes satis- 
fying Hq < k < eH(tf)a(tf), where e is a dimensionless 
parameter smaller than one. This guaranties that those modes 
are well outside the horizon at the end of inflation. Choos- 
ing e < a , the term -r((Vx) 2 ) turns out to be sub-leading. 
Moreover, if the moduli field is light (i.e., a <C 1) — which is 
the case of interest, at the end of inflation — Eq. (0 becomes 
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and has as solution \k — Bka, where Bk is a constant. This 
means that \k is constant at the end of inflation. Thus, the 
term (x 2 ) is also sub-leading. As a consequence, in this situ- 



ation one has p x 
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At this point, in order to calculate the quantity (x 2 ) dur- 
ing inflation, we first consider the de Sitter phase, in which 
case if we only take into account modes that leave the Hubble 
horizon, we get [6] 
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and taking now the derivative which respect to the cosmic 
time, we obtain 
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This can be generalized in terms of the stochastic equation [4] 
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which solution, for massive chaotic inflation, can be explicitly 
calculated: 



<x 2 > 



2v _ ZH 2 «{t){Hv 



H 



4-2a 



(*)) 



8(2 - a)ir 2 m 2 



(13) 



Here we have assumed that, at the beginning of inflation, x is 
not a condensate (i.e., (x 2 ) = 0, initially). 



For light particles, we have a <C 1 and, at the end of infla- 
tion, we get 
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which is much smaller than the result obtained in the pure de 
Sitter case 0] 
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because m x <C m. 



3. The moduli problem arising from gravitational particle 
production. — As we have seen, the gravitational particle pro- 
duction rate in the inflationary epoch could be very high, al- 
lowing for the possibility of creating a condensate (moduli 
field). Since at the end of inflation this field behaves as non- 
relativistic matter, it eventually dominates the energy content 
of the universe, thus violating standard nucleosynthesis pre- 
dictions, what will happen unless the amplitude of the mod- 
uli field is sufficiently small. The most stringent bound for 
successful nucleosynthesis prediction requires — provided the 
moduli field never dominates the expansion of the universe — 
that the ratio of the number density of light moduli produced 
during inflation, namely n x , to the entropy of the universe 
afterreheating, namely s, must satisfy the following inequal- 
ity 
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for moduli masses of the order of 10 2 — 10 3 GeV. 
In Ref. [3], the following expression was obtained: 
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where Tr is the reheating temperature, and Hq the value of the 
Hubble parameter at the beginning of inflation. (This formula 
comes from Eqs. (7) and (33) of Ref. [4]). On the other hand, 
in Ref. [4], for light moduli fields, the result was: 
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where t R is the time at which the reheating ends. 

If we consider the following realistic values: m 
(which ensures that the model does not overproduce aniso- 
tropics in the CMB signal), H — 1CT 4 A/ P (to obtain 10 4 
e-folds), and m x ~ 10~ 16 Af p (since we assume that the light 
moduli has a mass in the range of the gravitino one), then 
Eq. $T7\ yields the constraint Tr < 10~ 4 GeV, which is in- 
deed an abnormally small reheating temperature. On the other 
hand, using the relation 3H 2 (tR)M 2 = 71 J Q R , where g is the 
number of light degrees of freedom (for instance, g — 106.75 
for the standard model), Eq. $M yields T R > 10 16 GeV, 
which is in contradiction with the strong requirement that the 
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grand unified theory symmetries are not restored — in order to 
prevent a second stage of monopole production (see, for ex- 
ample, Hob- 

With the aim to understand these results we will now review 
the derivation of Eqs. (V7[ and (TTsT l, To obtain the first one, 
the authors assume that reheating occurs after the beginning of 
the oscillations of the moduli field. This gives the following 
bound to the reheating temperature: Tr < yJm x M p ~ 10 9 
GeV, which is needed to avoid the gravitino problem Isl fioll . 
They also assume that the field \ does not oscillate and that 
it (almost) does not change its magnitude up to the time, t, 
when H 2 (t) becomes smaller than m 2 /3. There, the ratio 
of the energy density of the moduli field with respect to the 
energy density of the infiaton field, at time i, is 
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where (x 2 ) is calculated at the end of inflation, since it is as- 
sumed that the amplitude of the moduli field almost does not 
change along the time interval (tf,t). 

From the end of inflation to the reheating time our patch 
of universe evolves as if it were matter dominated and, thus, 
the ratio (fl3T l is constant up to the reheating time. Then, 
as a consequence of the facts that at reheating time one has 
P4>{tn) = — Jo > that the entropy of the produced particles is 

27T 2 oT^ 

s = — an d that the number density of light moduli pro- 
duced during inflation is given by n x ~ the expression 
(7) in Ref. [3] was obtained, namely 
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To get the value of (x 2 ) one can assume that the moduli field 
is nearly massless, and then use Eq. (fl4] >. which yields 
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Finally, inserting it into Eq. (f20]l, one gets Eq. (fTTI l. 

Within a different approach, in Ref. J3l immediate thermal- 
ization was assumed for the infiaton and for the moduli fields 
after inflation, what leads to 
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The value of (x 2 ) for a light field at the end of reheating is 
taken in [4] from Eq. (Tl4"l i. what leads to the final result 
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4. Universal solution to the moduli problem. — After careful 
investigation, we here identify a problem which roots are com- 
mon to both situations — even if, in principle, the approaches 



look so different, even more the results obtained in the two 
cases. Indeed, in Ref. ||3J it is assumed that the amplitude of 
the moduli field does not change from the end stage of infla- 
tion to the beginning of the oscillations. Similarly, in Ref. fl] 
the assumption is done that the solution ( Tl4l > still holds after 
inflation. What we will here defend is that, immediately after 
inflation (H(tf) ~ m), when the universe behaves as if it 
was matter dominated, the moduli actually evolve like non- 
relativistic matter. Then, for t £ one has 
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Note that this expression is essentially the same as Eq. (3.3) 
of |8|], the crucial difference with Eq. ( l20b being the factor a 
in the numerator. To evaluate (x 2 ) at the end of inflation we 
use Eq. ( fT4b , and obtain 
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which, for realistic values: m ~ 10 e M p , Ho ~ 10 4 M p , 
and m x ~ 10 _16 M p , yields the perfectly consistent result 
Tr < 10 16 GeV. 

The main reason which leads us to assume that, immedi- 
ately after inflation has ended, the moduli evolve like non- 
relativistic matter is that the solution to Eq. ([6j, in the massless 
case, in a matter-dominated universe, is given by 111 llll2ll 
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Thus, for the infrared modes, one has 
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(22) which means that (x 2 ) 



As a consequence, for light 



moduli it is most reasonable to assume the same kind of beha- 
vior in the matter-dominated phase. 

Another way to understand this very important point, 
namely the fact that the moduli must evolve like non- 
relativistic matter, comes actually from the very requirements 
of massive chaotic inflation: the moduli field satisfies the con- 
servation equation 
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being p x the pressure. Therefore, as it occurs that in massive 
chaotic inflation the universe enters, at the end of the inflation 
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epoch, into a matter-dominated phase, one has that H(t) 



and thus the solution to Eq. (130b is p x 
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and p x ~ 0. 



Note, as a consequence, that although the mass of the field 
is rather small, it just cannot evolve like relativistic matter 
since it must necessarily satisfy Eq. (l30l l. while it turns out 
that relativistic matter will only satisfy this equation when 
H{t) ~ tjj (the scale factor would need to go as a(t) ~ t 1 ! 2 , 
and this is clearly not the case). 

To finish, a couple of remarks are in order. 

1 . If at the beginning of inflation one takes the energy dens- 
ity to be of order M p — which is the initial value considered, 
e.g., in (HQ— one gets H ~ M p and Eq. ([27) yields the 
constraint Tr < 1 GeV, which is indeed a very small value for 
the reheating temperature. Thus, in order to obtain a coherent 
result, one is necessarily compelled to assume that, initially, 
the Hubble parameter is smaller than M p . 

2. Numerical studies have been carried out with the purpose 
to calculate gravitational particle production. For instance, in 
Ref. IU5I1 a diagonalization method was used to calculate the 
production of minimally coupled light particles. As explained 
by the authors of that paper, such calculation is technically 
very difficult, for very small masses. Actually, what we do 
believe it really happens is that the method employed there 
only works well for massive conformally coupled particles 
Ill6lll7ll . In fact, the Bogoliubov beta coefficient is given by 



1 



\ Wx'Mt 



■ ul(v)\xk(v)\ 



1 



(31) 



where 



Xk{v) = a (v)Xk(v), 
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Note that, in the minimally coupled case £ = 0, one has 

ul(ri) = k 2 + a 2 {rj)(m 2 - 2H 2 - H). 
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Then, if the moduli is very light, for infrared modes this fre- 
quency is not well defined. Actually, in the de Sitter phase and 
for modes well outside the Hubble horizon, one obtains 



wgfa) = k 2 + a 2 (n)(m 2 -2H 2 ) 
~ k 2 - 2a 2 (r])H 2 < 0. 
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And, since these modes turn out to be the ones which are most 
important in order to calculate the relevant quantities [6, 18], 



we would finally conclude that the diagonalization method is 
not the most appropriate one to be used in this case. What 
is more, in order to compute the desired quantity, the authors 
calculate p x /m x numerically. However, p x is known to be a 
divergent quantity, if the field is not conformally coupled (see, 
e-g-, 11911 ). It must necessarily be renormalized previously, 
what renders utmost difficult to develop a rigorous procedure 
in order to perform such calculation numerically. 
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